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Abstract 

The evolution of a soliton star filled with fermions is studied in the 
framework of general relativity. Such a system can be described by 
the surface tension a, the bag constant B, and the fermion number 
density po . Usually one of these parameters prevails in the system and 
thus affects the spacetime inside the soliton. Whether it is described 
by Friedman or de Sitter metric depends on the prevailing parameter. 
The whole spacetime is devided by the surface of the soliton into the 
false vacuum region inside the soliton and the true vacuum region 
outside, the latter being described by the Schwarzschild line element. 
The aim of this paper is to study the equations of motion of the domain 
wall in two cases. In the first case the de Sitter metric describes the 
interior in the first case, and in the second case it is replaced by the 
Friedman metric. In both of them the Schwarzschild metric is outside 
the soliton. From the analysis of obtained equations one can draw 
conclusions concerning further evolution of a soliton star. 



In this paper we study the evolution of a soliton stellar configuration. Let 
us consider the SU C colour QCD plus a phenomenological colourless scalar 
field $ describing the colour confinement. The whole system is described 
by the scalar field $, the fermion quark field ipf, the gluon field A ^ and the 
gravitational field g^ u . The generalized QCD Lagrangian function takes the 
form: 

C = -\^~tl)F^F^ + \d^d^-U{^) (1) 

V 

and consists of two parts; the basic QCD Lagrangian supplemented by the 
Lagrangian of the dynamical scalar field. In the Lagrangian function (0) R 
is the curvature scalar and 

K = 8llG ~ "TTT~ (2) 
1V1 Pl 

G the gravitational constant and {/($) the potential. Due to nonlinearity 
of the potential U(<f>) the Lagrangian (||) leads to nontopological soliton 
solutions. The potential possesses the polynomial form of the fourth order 

i 

£/($) = |$ 2 + + + 5. (3) 

The constants a, 6, c are parameters which are adopted in the restriction to 
enable fitting of the static properties of hadrons. The bag constant B is 
a measure of the pressure of the physical vacuum on the perturbative one. 
The potential possesses two minima, one at $ = which is a local minimum 
associated with the perturbative vacuum state and the second, the absolute 
minimum, at 

, 36 M _ Sac.u , A . 

corresponding to the physical vacuum. The form of the potential is presented 
on Fig.l. The soliton solution exists between two limiting cases B = which 
represents the degenerate vacuum state and B = B max when the local mini- 
mum and maximum overlap and an inflection point appears. This is the mo- 
ment when the perturbative and physical vacuum states become metastable 
and also the limiting case in which the soliton solution still exists. 
The action 

S = J d A x^jC (5) 
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with C defined as (^) leads to the appropriate equations of motion. We are 
looking for spherically symmetric solution with the metric 

ds 2 = g^dxPda? = -e u dt 2 + e x dr 2 + r 2 [dB 2 + sin 2 6d^ 2 ) (6) 

Variation of the action (|5|) with respect to metric g^ u gives the four-dimensional 
Einstain equations 

Rfiu — -^9^R = (7) 
with the energy-momentum tensor T^ v 

T^ = T^)+T^) (8) 

where 

T^($) = d^d,,® - gta ,C{$) (9) 

is the scalar field energy- momentum tensor and T^ip) its fermion counter- 
part. 

The soliton bag model in the thin-wall approximation describes the quark 
star model with a surface tension. The whole spacetime is devided by the 
surface of the soliton into the false vacuum region inside the soliton and the 
true vacuum region outside, the latter being described by the Schwarzschild 
line element. The aim of this paper is to study the dynamics of the domain 
wall in two cases. In both of them the Schwarzschild metric is outside the 
soliton. 

Subscripts ± refer to values of the corresponding coordinates in the region 
outside and inside the soliton. 

ds 2 + = -(1 - —)dt 2 + + dr 2 + r 2 (dd 2 + sin 2 9d0 2 ) (10) 

r v r+ > 

where r g = 2GM with M being the total soliton mass. 

The de Sitter metric describes the interior of the soliton in the case when the 

value of the cosmological constant A dominates. 

ds 2 = -(1 - -Ar 2 )dt 2 + \ — ^-dr 2 + r 2 (d6 2 + sin 2 Bd$ 2 ) (11) 

v 3 ' (1 - |Ar 2 ) v ' v ' 

In the second case when the fermion number density p is t ne prevailing 
parameter, the interior metric is replaced by the Friedman solution M. The 
line element has now the form 

dv^ 

dl 2 = + r 2 [dd 2 + sin 2 6d(j) 2 ) (12) 
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where a is a curvature radius. After performing the following substitution 
r = a sin x where x changes from to n, the line element takes the form 

dl 2 = a 2 {d X 2 + sin 2 X (d9 2 + sin 2 Odcj) 2 )} (13) 

For space with constant, positive curvature and for a = ar g where a is the 
dimensionless parameter, the metric can be written in the following way 

ds 2 = -dt 2 + a 2 r 2 d X 2 + a 2 r 2 (d6 2 + sin 2 6d<f) 2 ). (14) 

Because 

dv 

Tgdx = —^2 (15) 
1 r 2 



finaly one can get the metric 



ds 2 = -dt 2 _ + a 2 



' dr 2 _ 



r 2 
r 2 



+ r 2 _(d9 2 + sin 2 9d(f) 2 ) \. (16) 



The time evolution of the scale factor a inside the soliton star resembles cos- 
mological evolution of the universe but on a different scale. Inside the soli- 
ton star we have a nonvanishing cosmological constant A = kB and fermions 
which at rough approximation are treated as dust density p Q . The space-time 
in the soliton star will evolve according to the Einstein equation which gives, 
as in cosmology, the scale evolution equation 

4 1 
a 2 + 1 = -TvGp a 2 + -Aa 2 (17) 

o o 

It is possible to interpret the time evolution of the metric inside the soliton 
as a movement of the fictitious particle with the vanishing energy under the 
influence of some effective potential 

l -a 2 + U eff {a)=E (18) 

For E = 

Ueff(a) = -- -g-(Po^)- - -g-£a 2 (19) 

where the dust mass M d = p V with V = 2n 2 a 3 . The "steady state" solution 
is obtained for 

GM d i M d i 
ao = ( ^r )3 = (20) 
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It is convenient to introduce the dimensionless potential 

1 . A 1 



2-^+2^ 



U eff(x) 

which depends only on one parameter A and a 

37ran 



a x 



(21) 



(22) 



The critical point A c — | represents the maximum of the potential U e ff (x)Fig.2 
presents the form of the potential function U e ff(x) for two values of parame- 
ters A; A = A c and A — |. Results obtained for different values of parameter 
B are presented in Table 1. 



B = (O.lGeV) 4 


B = (IGeVf 


ao = 1036 km 


a = 0.48 km 


A = 0.001524 


A = 0.3284 



The metric evolution is determined by the critical value of parameter A c 
which produces the critical scale a c . If A > A c , the cosmological constant A 
can be neglected and the evolution is exactly the same as in the Friedman 
universe with gravitational collaps as a result. For the soliton star with 
B = 56MeV/ fm 3 we obtain p c ~ 10 11 g/cm 3 which is below the neutron star 
density. If A < A c , the cosmological constant A dominates and space-time 
expands. The empty soliton will expand exponentially according to the de 
Sitter solution with the Hubble constant 



H 



I8tiG 



B. 



(23) 



The bag constant B=56 MeV//m 3 determines the value of the parameter tn 



- ~ 10> S 



lday. 



(24) 



In order to derive the equations of motion for the surface of a soliton star we 
have to obtain the junction conditions between the inner and outer region 
of the soliton. Let E be a three-dimensional spacetime hypersurface swept 
out by the soliton surface. It is convenient to choose the Gaussian normal 
coordinates in the neighbourhood of S. Coordinates of any point p belonging 
to the neighbourhood N of E are given by x M = (xj, rj), Xi = (r, 9, ip) where ip 
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and 6 are angular variables, r the proper time of an observer comoving with 
the domain wall, 77 is the distance from X to p along the geodesies which 
are orthogonal to E ||. We also define as a unit vector field normal 
to 77 = const hypersurface. The three dimensional metric intrinsic to the 
hypersurface is 

hfj,u = 9iiv ~ £,fj,£i> (25) 

where g^ v is the four metric of the spacetime. According to Israel 0, in this 
case the 3-tensor 7^ defined as 

7ij Kj'j A / (26) 

where Ky is the extrinsic curvature 3-tensor of r] = const hypersurface. This 
tensor can be calculated using the following relations 



- -ns - It- < 27 > 



We are interested in the situation when the energy-momentum tensor of 
the four-dimensional spacetime has a 5-function singularity on the hypersur- 
face and can be expressed in the following form |5|] 

Tpv = S^{ Xi )5{rj) (28) 

with S^v as the surface energy tensor. Taking this into account the extrin- 
sic curvature possesses a jump discontinuity across the hypersurface. 
Adopting the Gauss-Codazzi formalism []|] the 7^' tensor takes the form 

7i . = -8ttG0% - l - 9ij S) (29) 

where is defined as the integral of the energy-momentum tensor through 
the thickness of the wall. 



= lim / Tijdx 1 (30) 


We will limit our consideration to the following form of the surface energy 
tensor || 

Sp, = O^U^Uv - C(t)(V/ + u v u v) ( 31 ) 
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where a is the surface energy density of the domain wall, £ the surface tension 
and w M the four-velocity of the wall. For the domain wall o = ( and 

= -ah^ (32) 

In order to describe the dynamics of the wall which separates the interior 
and the exterior of the soliton star we have to calculate the appropriate 
components of the extrinsic curvature; K~t - the extrinsic curvature for the 
Schwarzschild metric and K^j for the metric which describes the space inside 
the soliton. Taking into account the form of S^ u ( |32| ) 

7 j = -AnGaS] (33) 

For the system described by the Lagrange function (Q) it is possible to inter- 
pret the Einstein equation having the following form inside the soliton 

Rpu ~ -^QnuR = ^Bg^ + kT^i!)) (34) 

or 

Rfiu - ^9^R - A 9vu = kT^i/)) (35) 

with A = kB as the Einstein equation with the cosmological constant differ- 
ent from zero and with the fermion energy-momentum tensor. In the case 
when T^ixp) = (empty bag) the Einstein equation is reduced to the form 
with the solution expanding according to the de Sitter law. It is interest- 
ing to present the dynamical equation of the soliton surface. Owing to the 
spherical symmetry of £ the metric on the domain wall can be presented in 
the following way 

d s 2 = -dr 2 + R(T) 2 (d6 2 + sin 2 Odcf) (36) 

where r is a proper time along the world-line (6, 4> = const) of a dust particle. 
We adopt this metric as a matching condition of the interior and the exterior 
metric on the soliton surface. Now using the method presented in Israel's 
paper || we obtain the dynamical equation for the soliton surface 

(1-l-AR 2 + R 2 )^ -(l-^ + R 2 ) =A<kcjR (37) 
3 R 

In all these equations R = 4p. This equation agrees with that obtained by 
A.Aurilia 0. Once more the surface dynamics can be illustrated as fictitious 
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particle motion with the total energy equal to zero in the effective potential 
field. 

\k 2 = -V eff (R) (38) 
where V e ff(R) can be expressed as 

Veff(R) = 1 -{l-R 2 (^ + 2irGaf 
M B M 2 

~R [ ^G^ + l) ~ IGttV^ 1 m 

substituting R = xr g , r g being the gravitational radius, the effective potential 
V e ff(R) takes the following form 

Veff = ^{(l-x 2 r 2 (^- + 2nGa) 2 
IB 1 

2^ ( 6^2 + X) " (8%Gr g( x) 2 x^ (40) 
We can also present the mass M of the soliton 

M = \nR*B + ^R 2 a{\ - \tiGBR 2 + i? 2 )^ - 8ir 2 Ga 2 R s (41) 

The subsequent terms in the expression for the mass are the volume energy 
and the surface energy with the relativistic correction and the surface-surface 
binding energy. Fig. 3 presents the form of the potential V e ff(x) which de- 
pends only on one variable x. The soliton described in this paper is not 
stable and such a configuration tends to obtain the equlibrium state by the 
emission of particles. Phenomenological description of this situation can be 
made by introducing the term identified as friction. 

In the vicinity of point x c which corresponds to the maximum of the potential 
V e ff(x), the function V e ff(x) can be approximated in the following way. 

1 1 

On Fig. 3 the function V e ff ta (x) is represented by the broken line. Substituting 
y = x — Xo one can write 

y 2 = u 2 {y + a 2 ) (43) 
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where a 2 = The general solutions of this equation have the following 
form 

y x = i{e-'°e^(-a 2 + e 2ta e - 2uJt )} (44) 

y 2 = i{e-'°e^(-a 2 + e 2t °e 2ujt )} (45) 
where y\ represents collaps and y% expansion of the soliton. For t = 

Vo = \{e- t0 {~a 2 + e 2t °)} = \{~a 2 e^ + e to } (46) 
so for this moment 

x p = x + ^{-a 2 e~ t0 +e t0 } (47) 

Substituting e'° = az 

x p = x + ~(z 2 -l) (48) 

The mass M of the soliton is determined by the initial value x p . The dimen- 
sionless parameter m = 

4 

m = -Tibx^ + 4-nsxl (49) 
3 

The equation of motion takes the following form 

x + bx = -^tl* = uj 2 (x - x c ) (50) 
ox 

We can analyze the solution of the equation floTf ) in momentum space Fig.4. 
Solutions where p > 0, p = x, describe the trajectories representing the ex- 
pantion of the whole system whereas for solutions with p < the trajectories 
lead to the collaps. Thus the values of the initial conditions are responsible 
either for the expansion or the collapse of the soliton. There is only one 
critical trajectory for which the whole system, after an infinite time, tends 
to the stable configuration M = O.92M and R = 16.33 km. Different types 
of trajectories are presented on Fig.4. The vertical stright line marks the 
position of the Schwarzschild horizon, two hyperbolas correspond to the tra- 
jectories if there is no friction in the equation of motion, broken line presents 
the trajectory which passes very close to the critical one. The stright line 
shows the position of the critical trajectory. The analysis of this figure leads 
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to the conclusion that the system reaches stable configuration before reaching 
the Schwarzschild horizon. Such a configuration corresponds to the situation 
when the soliton obtaines the state of equilibrium between the emission of 
particles from the surface and the de Sitter expansion. Fig. 5 shows the nu- 
merical solution of the equation of motion for the potential function V e ff(x) 
and this solution confirms the analysis which was carried out on the basis of 
Fig.4. 

Let us now concentrate on the situation when the cosmological constant A 
is neglected inside the soliton. It has been already mentioned that this very 
solution corresponds with the case when the bag constant B = and rep- 
resents the degenerate vacuum state. The confinement is only due to the 
surface tension a. This situation appears in the case of the Lee-Pang star 
0. The necessary and sufficient conditions for reality of this system is both 
the existence of the nontopological soliton which is guaranteed by the scalar 
field the fermion field and the gravitational field. The soliton star in this 
case consists of a large interior and a much thinner surface shell of the width 
~ ^ — 1), where \i is the mass of the scalar field $. For the value of surface 
tension a = ^(30GeV) 3 the following parameters of the soliton star was ob- 
tained: massM ~ 1O 1 3M and radius R ~ 1 light year. 
The exterior geometry is still described by the Schwarzschild line element 
whereas the interior geometry is described by the Robert son- Walker metric. 
If fermions inside the soliton are treated in the rough approximation as dust 
with the density po> the energy momentum tensor T^ u takes the form 

Tpv = Pqu^Uv (51) 

where is the four-velocity vector of the particle. Once more the Israel's 
equation for the motion of the soliton surface separating the two metrics was 
obtained 

h-\l + h&)* -(l-?^- + R 2 )s =AiraGR (52) 

R 

2 

where h = , a is the scale factor. In order to study the dynamics 

of the soliton surface we identified the motion of the domain wall with the 
motion of a unit mass particle with vanishing total energy moving under 
the influence of the potential. The form of the potential function W{$) is 
presented on Fig.6. Closer analysis of this function reveals that once again if 
we introduce the friction the whole system behaves just like in the de Sitter 
case. One can notice the similarity with the results obtained in the previous 
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paragraph concerning the analysis of the function V e ff(x). However, this 
time the critical trajectory reaches the Schwarzschild horizon R = r g . Thus 
in this case the soliton collapses to the black hole. 

1 Conclusions 

The analysis of the soliton surface dynamics has been performed in this pa- 
per. It may be concluded that the star either collapses into a black the whole 
system is characterized by the following parameters the bag constant B, the 
surface tension a, and the fermion number density po- If the cosmological con- 
stant A is the prevailing parameter the space time inside the soliton can be 
described by the de Sitter solution. When p prevails the Friedman solution 
is obtained. It turned out in the de Sitter case there exists the stable con- 
figuration described by the macroscopic parameters: the mass M = 0.92M Q 
and the radius R = 16.33 km. The radius of this configuration is bigger than 
the gravitational one and chartacteristic of the whole system resambles the 
ones that can be observed in neutron stars. In the case of Friedman solution 
the soliton which is characterized by the mass M ~ 1O 13 M collapses exactly 
at the Schwarzschild horizon, so in the result we obtain the black hole. 
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Figure captions 



Fig.l. The form of the potential function 

Fig. 2. The effective potential U e ff(x) for two values of parameter A: the 
continues line corresponds to the A = A c wherease dotted line to A — |. 
Fig. 3. The form of the effective potential V e ff(x) - continues line, dotted 
line indicates the possition of the function V e ff ta (x). 

Fig. 4. The trajectories for different types of solutions presented in the phase 
space. 



Fig. 5. Numerical solution of the equation of motion |50 
Fig.6. The effective potential W(§). 
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Figure 2: The effective potential U e jj(x) for two values of parameter i 
continues line corresponds to the A = A c wherease dotted line to A = 




Figure 3: The form of the effective potential V e ff(x) - continues line, c 
line indicates the possition of the function V e ff >a (x). 
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Figure 4: The trajectories for different types of solutions presented in the 

phase space. 
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